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Abstract. Dirichlet averages of multivariate functions are employed 
for a derivation of basic recurrence formulas for the moments of multivariate 
Dirichlet splines. An algorithm for computing the moments of multivariate 
simplex splines is presented. Applications to hypergeometric functions of 
several variables are discussed. 

Introduction. 

In [9], H.B. Curry and I.J. Schoenberg have pointed out that univari- 
ate S-splines can be constructed from volumes of slices of convex polyhedra. 
An extension of this idea to the case of multivariate splines is due to C. de 
Boor [1]. Since the geometric construction is too complicated to be used in 
numerical computations, some basic recurrence formulas for these functions 
have been found (see [3,8,10,11,14,17,18,19,22,23]). Multivariate S-splines 
(also called simplex splines) have been studied extensively over the past 
thirteen years by many researchers. These functions have been found useful 
for some applications of data fitting, computer aided geometric design, and 
mathematical statistics. In [29] the author addressed some new problems 
where the simplex splines could play a prominent role. 

Further generalizations of multivariate simplex splines appear in [19] 
and [13]. Here and thereafter we call these functions Dirichlet splines. Our 
choice of terminology is motivated by the fact that the distributional defini- 
tion of the Dirichlet splines (see (2.2)) involves the Dirichlet density function. 
It has been demonstrated that this class of splines is well designed for some 
problems of mathematical statistics (see [13], [19], and references therein). 
Applications to the theory of multivariate convex functions are reported in 
[26, 27]. 

Recently we have noticed that there is a simple relationship between 
univariate Dirichlet splines and some special functions such as 
-R-hypergeometric functions and confluent hypergeometric functions 



(see [28]). In this paper, we present some results for multivariate Dirichlet 
splines together with their applications to special functions of several vari- 
ables. 

The outline of this paper is as follows. In Section 2 we introduce 
notation and definitions which will be used throughout the sequel. In the 
third section we give a definition and basic properties of a Dirichlet average of 
a multivariate function. Also, we give a recurrence relation for these averages. 
In Section 4 we give two results which play a crucial role in our subsequent 
considerations. Moments of the class of splines under discussion together 
with two moment generating functions are presented in Section 5. Therein, 
we also give an algorithm for computing the moments of multivariate simplex 
splines. Applications to hypergeometric functions, including Appell's F 4 and 
Lauricella's Fb, are discussed in Section 6. In the same section, we give a 
recurrence formula, two generating functions, and an inequality for Lauricella 
polynomials. 

2. Notation and Definitions. 

Let us introduce some notation and definitions which will be used 
throughout the sequel. By x, y, . . . , we denote elements of Euclidean space 
]R S (s > 1), i.e., x = (xi, . . . ,x s ). Superscripts are used to number vectors. 

s 

The inner product (or dot product) of x, y G M s is denoted by x-y = Yl x kVk- 

k=l 

For a given set Icl s the symbols [X] and vol s (X) mean the convex hull of 
X and the s-dimensional Lebesgue measure, respectively. We use standard 

multi-index notation, i.e., for (3 G Z+, \/3\ = (3\ + h f3 8 , (3\ = (3\\ ■ ■ ■ f3 8 \, 

a < 13 means a\ < . . . , a s < j3 8 (a G Z+), x 13 = xf 1 x^ 2 ■ ■ ■ x^ s . For 
r G Z + , and /3 E Z s + with \/3\ = r, the multinomial coefficient (V) is defined 



in the usual way 



( 




By 



n 



E n = {t = (£1, . . . ,t n ) G R n : tj > 0, for all j, ^t 3 < 1} 



we denote the standard n-simplex. Let R> represent the set of all positive 
real numbers and let b = (6 , ■ ■ ■ ,b n ) E The Dirichlet density function 

on E n , denoted by 4>b, is given by 



For X = {x ,... ,x n } C M s (n > s > 1) with vol a ([X]) > 0, the 
multivariate Dirichlet spline M(-\b;X) is defined by requiring that 



holds for all / G Co(lR s ) - the space of all multivariate continuous functions 
on W with compact support (see [13, 19]). Here dx = dxi---dx s , dt = 



dt\ ■ • -dtn, Xt = Ux 1 . When bo = • • • = b n = 1, 4>bif) = n\, and the 

i=0 

corresponding spline becomes a simplex spline. The latter spline will be 
denoted by M(-\X). When s = 1, we will write m(-\ b; Z) instead of M(-\b;X), 
where Z = {zq, . . . , z n } C K. In this case (2.2) becomes 



n 



Mt)=r(c)H[T(b t )]-H^-\ 



where t G E n , to = 1 — t\ — ■ ■ ■ t n , and c = bo + ■ ■ ■ + b 




n 




(fteCo(R)). 



3. Dirichlet Averages. 

The purpose of this section is two-fold. We give a definition of Dirich- 
let averages of multivariate functions. Next we prove a recurrence formula for 
these averages. This result has an immediate application in Section 5. For 
the reader's convenience, let us recall a definition of the Dirichlet average of 
a univariate function h G Co(K). Assume that the set Z = {z , . . . , z n } C K 
is such that mm{zj : < j < n} < max{zj : < j < n}. For b G the 
Dirichlet average of h, denoted by H(b; Z), is given by 



We list below some elementary properties of H(b; Z). 

(i) H(bo, ... ,b n ;zo,... , z n ) is symmetric in indices 0, 1, . . . , n (see [6, 
Thm. 5.2-3]). 

(ii) A vanishing parameter bi can be omitted along with the corresponding 
variable Zi (see [6, (6.3-3)]). 

(iii) Equal variables can be replaced by a single variable if the correspond- 
ing parameters are replaced by their sum (see [6, Thm. 5.2-4]). 

We now introduce the Dirichlet average of / G Co(K s ). For X C M s 
with vol s ([X}) > and b G M> +1 (n > s > 1), the Dirichlet average of /, 
denoted by JF(6; X), is given by 




(see [6]). Comparison with (2.3) yields 





5 



where Xt and <pb have the same meaning as in Section 2. Comparison with 
(2.2) shows that 

F(b;X)= I f{x)M(x\b;X)dx. (3.4) 

It is clear that the properties (i)-(iii) are also valid for the average T . 
In particular, property (iii), when applied to M(-\b;X) yields 

M(-\b;X) =M(-\x°,...,x°,...,x n ,...,x n ) 
(b — times) (b n — times) 

provided the 6's are positive integers. (See also [13, 19].) The spline on the 
right hand side of (3.5) is a multivariate simplex spline with coalescent knots 
(see [14] for a detailed analysis of this class of splines). 

Before we state and prove the first result of this section, let us intro- 
duce more notation. By ej (0 < j < n), we denote the jth coordinate vector 
in R n+1 . For / G Cq(M s ) define a function fj as follows: 

fj(x) = D x i_J(x) (3.6) 

Here D y f denotes the directional derivative of / in the direction y e 1R S , i.e., 

We are now ready to prove the following. 

THEOREM 3.1 ([!}). Let X = {x°, ...,s"}cK s (n>s>l)be such that 
vo\ 8 ([X}) > 0. Further, let / G C^R 3 ) and let the vector b G M> +1 be such 
that bj > 1 for some < j < n. Then the following identity 

{c-l)T{b;X) = (c-l)F(b-ej;X) + Fj(b;X) (3.7) 

is valid. Here Tj denotes the Dirichlet average of the function fj. 



REMARK. The proof presented below bears no resemblance to what was 
done in [7, Thm. 3]. In this paper, the author has established (3.7) using 
generalized Euler-Poisson partial differential equations. 

Proof. In order to establish the identity (3.7) we employ the following one 

(c - l)H(b; Z) = (c- l)H(b - ey Z) + Hy{b; Z). (3.8) 

Here Hj stands for the Dirichlet average of the function hj (u) = (zj —u)h'(u), 
h G Cg(R s ). The relation (3.8) readily follows from (5.6-13) in [6]. Applica- 
tion of (3.2) to (3.8) yields 

(c — 1) / h(u)m(u\b; Z)du = (c — 1) / h{u)m{u\b — ey Z)du 

Jr Jr ^ g^j 

+ / (zj — u)h! (u)m(u\b; Z)du. 
Jr 

We will lift (3.9) to the case of multivariate functions. To this aim we shall 
employ the following formula 

/ h(u)m(u\b; Z)du = / h(X ■ x)M(x\b; X)dx, (3.10) 
Jr Jr s 

where now Z = {\-x°, . . . , \-x n }, A e M s \{0}, and h(X-x) is a ridge function 
(or plane wave). Since the proof of (3.10) is similar to that presented in [23], 
p. 496, we omit further details. Application of (3.10) to (3.9) yields 



(c-1) / h(X-x)M(x\b;X)dx= (c- 1) / h(X ■ x)M(x\b - ey, X)dx 
Jr s Jr s 

+ [ (X-x j - X-x)h'(X-x)M(x\b;X)d. 
Jm a 



(3.11) 

x. 



We appeal now to the denseness of ridge functions (these functions form a 
dense subset in Cq(1R s )) to conclude that the above identity is valid for any 
multivariate function / e Cq(1R s ), see [20]. Substituting h(X ■ x) = f(x) into 
(3.11), we obtain the assertion and the proof is completed. □ 
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COROLLARY 3.2. Along with the hypotheses of Theorem 3.1, assume that 
for some < i, j < n, 1 < k < s, that x\ ^ 0, x J k ^ 0, and bi > 1, bj > 1. 
Then 

(c - 1) [F(b - ej ;X) - F{b - ef,X)} + T 3 {b- X) - F t (b; X) = (3.12) 

and 

(c - 1) (xl - 4)^(6; X) = (c - 1) [4^(6 - e, ; X) - x{F{b - e,; A)] 

(3.13) 

+ x t k r j (b;X)-x 3 k F i (b;X). 

REMARK. (3.12) and (3.13) are generalizations of Exercise 5.9-6 in [6] with 
the latter being an extension of Zill's identity for .R-hypergeometric functions. 

Since the proof of (3.12) and (3.13) follows the lines introduced in [6, 
p. 305], we omit further details. 

4. Auxiliary Results. 

Our first result reads as follows. 

PROPOSITION 4.1. Let p(x) be an affine function on M s . Then 

n 

p(x)M(x\b; X) = J2 w l p(x l )M(x\b + a; X) (4.1) 

provided the splines M(x\b + ef,X), < i < n, are continuous at x G M s . 
Here Wi = bi/c, i = 0, 1, . . . , n. 

Proof. We need the following identity for the Dirichlet density function [6, 
(4.4-8)]: 

U<l> b (t) = w^ b+ei {t) (4.2) 
(t e E n ; < i < n). Since t H \-t n = l, (4.2) gives 

n 

fo(t) = ^Witb+ejjt)- 
i=0 



Multiplying both sides by f(Xt) and next integrating over E n , we obtain by 
virtue of (2.2) 

/ f(x)M(x\b;X)dx = [ f(x)J2wiM(x\b + e t ;X)dx. 

JR S Jr s i=Q 

Hence, (4.1) follows when p(x) = 1. To complete the proof we utilize (4.2) 
again. Multiplying both sides by x\ and next summing over i, we obtain 

n 
i=0 

Here (Xt)i denotes the Zth component of Xt. This leads to the following 
integral relation: 



/ f(x)xiM(x\b;X)dx = [ f(x) V Wi x\M{x\b + e*; X) 
JR a Jr s \-i=o 



dx, 



which proves (4.1) when p(x) = xi, (1 < I < s). □ 

Micchelli [22] gave a different proof of (4.1) for simplex splines. For 
this class of splines, identity (4.1) is called the "degree elevating formula". 
A special case of (4.1) appears in [14]. 

For our further aims, we recall a definition of the .R-hypergeometric 
function in the real case. Let the set Z = {zq, . . . , z n } C R n+1 be such that 
i [Z\. Further, let b G M> +1 . The .R-hypergeometric function R_ a (b;Z), 
(o 6 1), is given by 

R. a (b-Z)= [ (Zt)- a <p b (t)dt 

(see [6]). When —a e N, the restriction ^ [Z] can be dropped. Comparison 
with (3.1) shows that the R- a is the Dirichlet average of the power function 
u~ a . It is worthy to mention that the Gauss hypergeometric function 2-F1, 
Lauricella's hypergeometric function Fx>, the Gegenbauer polynomials, and 



the elliptic integrals in the Legendre form can all be represented in terms of 
the function R- a - Combining (4.3) and (2.3) gives 

R- a (b;Z)= / u~ a m(u\b; Z)du. 
Jr 

For future use, let us record a very useful formula for 
-R-hypergeometric functions (see [6], Thm. 6.8-3) 

n 

R_ a (b; Z) = Y[zJ bi R a - c (b; Z" 1 ), 

3=0 

where Z' 1 := {z^ 1 , . . . , z' 1 } (zj > 0, for all j; c ^ 0, -1, -2, . . . ). This 
important result is commonly referred to as Euler's transformation. 

We close this section with the following: 

PROPOSITION 4.2. Let a e R and let the vector A G I s be such that 
\-x* < 1, j = 0, ... ,n. Then 

/ (l-\-x)- a M(x\b;X)dx = R- a (b;Y), 

where 

Y = l-X-X = {(l-X- x°), . . . , (1 - A • x n )}. 

Proof. Substituting Z = Y into (4.4) we obtain 

R_ a (b;Y)= [ u~ a m(u\b; 1 - A • X)du 
Jr 

= [ (1 - u)- a m{u\b-\- X) 
Jr 

= [ (1- A- X)- a M(x\b;X)dx. 
Jr s 

In the last step we have used (3.10). □ 



When a = c, (4.6) becomes Watson's identity (see [31]) 

„ n 

/ (l-X-x)- c M(x\b;X)dx = T\(l-\-x j )- b i. 

Jr* j=0 

The above identity follows by applying (4.5) to the right side of (4.6) and 
using R = 1. 

An alternative proof of (4.8) appears in [13] (see also [19] for some comments 
concerning this identity). Subject: file 3 

5. Moments of Multivariate Dirichlet Splines. 

A motivation for the investigation of the moments of Dirichlet splines 
has its origin in two mathematical disciplines. It is well known that the spline 
M(-\b;X) is a probability density function on R s . We feel that the results 
of this section can be applied to some problems in mathematical statistics. 
A second area of possible applications is the theory of special functions. 
We have already mentioned that some important special functions can be 
represented by the .R-hypergeometric functions. For particular values of the 
a parameter and the b parameters in (4.4), this function becomes a complete 
symmetric function. For particular values of the ^-variables, (4.4) gives an 
integral formula for the (/-binomial coefficients (Gaussian polynomials). (See 
[25] for more details). 

In this section, we derive recurrence formulas for the moments of mul- 
tivariate Dirichlet splines. Also, we discuss implementation of these results 
in the case when b = (1, . . . , 1) e M n+1 . For related results when s = 1, see 
[24]. We employ the multi-index notation introduced in Section 2. 

For e M s , we define the moment of order \/3\ (\/3\ = /?i H + f3 s ) 

of M(-\b;X) as follows: 

m p {b;X)= [ x l3 M{x\b]X)dx 
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provided S ^ [X], S - the origin in M s . When (3 G Z^_, this restriction is 
nonessential. In the case of the simplex spline, we shall omit the vector b and 
write mp(X) instead of mp{b;X). Also, let di stand for the Zth coordinate 
vector in M s . 

We are now ready to state and prove the following. 

THEOREM 5.1. Let the weights wo, . . . ,w n be the same as in Proposition 
4.1. Then 

n 

mf3(b;X) = ^2w i m f3 (b + e i ;X), (5.2) 

i=0 
n 

mf3 +dl {b; X) = Wix\mp(b + a\ X), (5.3) 

i=0 

for all / = 1,2,... , s. Moreover, if vol s ([X]) > and bj > 1, for some 
< j < n, then 



(c + |/?| - l)m p (b- X) = {c- l)m p (b -e ;X) + J2 Pi4 m P-di (&; X). (5.4) 

i=i 

If for some < i, j < n, 1 < k < s, x\ ^ 0, x 3 k ^ 0, and bi > 1, bj > 1, then 

s 

(c - 1)^(6 - e i5 X) -m p (b- e,; X)] + ^ /3 fc (4 - 4)™/3-d* (6; X) = (5.5) 

fe=i 

and 

(c+ - 1)(4 - xfompfa X) = (c - l)[4™/?(& - ey, X) - x{m p {b - a; X)} 

s / r f>\ 



i=i 



i T 3 
k X k 

I ry.3 



x l x x 



where W^i = det 

REMARK. When bj = 1, formula (5.4) holds true provided n > s. 

Proof. In order to establish the recursion (5.2) and (5.3), we substitute 
p(x) = 1 and p(x) = x\ respectively, into (4.1) and next integrate over R s . 
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For the proof of (5.4), we utilize formula (3.7) with f(x) = x@. The resulting 
equation, together with (3.6) and (5.1), yields the assertion. Formulas (5.5) 
and (5.6) follow immediately from Corollary 3.2 with f{x) = x&. □ 

We now give two moment generating functions. The first generating 
function involves the confluent hypergeometric function S. Following [6, 
(5.8-1)], we define 

S(b;Z)= exp(Zt)<f> b (t) dt, (5.7) 
b e K> +1 , Z = {z , . . • , z n }. Use of (3.1) and (3.2) gives 



S(b; Z) = / exp(u)m(u\b; Z) du. 
Jr 



Letting Z = A • X = {A • x°, . . . , A • x n } : A e M s \{0}, and next using (3.10), 
we arrive at 

S(b;X-X)= I exp(X-x)M(x\b;X)dx. (5.8) 
Jr s 

To obtain the first moment generating function, we expand exp(X ■ x) 
into a power series. Applying the multinomial theorem to powers of A • x 
and next integrating the corresponding power series one term at a time, we 
obtain by virtue of (5.8) and (5.1) 

S(b;X-X) = J2jrn J (b;X) (5.9) 

where the summation extends over all multi-indices j G 

It is worthy to mention that the hypergeometric function S(b; X ■ X) 
can be expressed as a divided difference of exp(z) provided that 
bo, . . . , b n e Z + . We have, 

kl S(b; X ■ X) = [A • x°(b ), . . . , A • x n (b n )]e z (5.10) 
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where k = b H b n — 1. Here the symbol A • means that the knot set 

A • x 1 is repeated 6j times. Formula (5.10) readily follows from (5.8), (5.7), 
and the Hermite-Genocchi formula for divided differences. 

A second generating function is given by 

R. a {b- Y) = J2 A J rnj(b; X) (5.11) 

(| A • x l \ < 1, for all i), where the summation extends over all multi-indices 
j G Z s + . Here, a G M, the set Y is given in (4.7), and (a,/) stands for the 
Appell symbol, i.e., (a, 0) = 1, (a, /) = a(a + 1) • • • (a + 1 — 1), I e N. In order 
to establish (5.11) we expand (1 — A-x) _a into a power series and next utilize 
the multinomial theorem to obtain 

i=o \j\=i J - 

To complete the proof we substitute (5.12) into (4.6) and next inte- 
grate term by term. Applications of (5.9) and (5.11) are discussed in the 
next section. 

We shall now turn our attention to the case of multivariate simplex 
splines. To this end, let f3 G Z^_. In the case under discussion, the formulas 
of (5.3) and (5.4) of Theorem 5.1 take the form 

1 

m p+dl (X) = Y. x \ m ^ Xi ) ( 5 - 13 ) 

Tl T 1 . 

1=0 

s 

(n+\p\)m /3 (X)=n[m (3 (X j )}+J2fo xi i m P-di( X )' ( 5 - 14 ) 

l=i 

1 < j < n; 1 < / < s. Here X { = X U {x 1 }, and X, = X\{x { }, < i < n. 
The set X 1 appears on the right hand side of (5.13) because of (3.5). A closer 
look at (5.14) shows that the recursion is in two directions. That is, given 
X = {x°, . . . , x n } C 1R S (n > s), to compute mp(X), we need the moment 



of order \(3\ for the knot set consisting of one less vector than X, and also s 
moments of order \(3\ — 1. 

Define a set X k = ... , x k }, k = s, s + 1, . . . , n, and note 
mp(Xk) = 1 when |/?| = 0, k = s, s + 1, . . . , n. To employ (5.14), we must 
precompute certain moments of the form mp(X s ), \/3\ > 0, and m,;, (X^), 
k = s, . . . , n. Let us note that 

_ 1 k 

m dl {X k ) = 1 —^J2 x i J 

3=0 

follows immediately from the defining equation (2.2). 

In order to compute the moments m i a(X s ), \(3\ > 0, we first introduce 
some new notation and next appeal to the proposition that follows. 

Let t = (to,t!,... ,t n ) G R n+1 with ,t n ) G E n , and t = 

n 

1 — ti. Also, let / = (Z , . . • , l n ) G Z" +1 , with m = |/|. Then we define 

i=l 

the Bernstein polynomial by 

For given coefficients {pi}, \l\ = m, we shall call any polynomial of 
the form 

q(t) = Pl B i m (t) 
\l\=m 

a Bezier polynomial. It is well known that any such q may be stably and 
efficiently evaluated using deCasteljau's algorithm (see [2], [12]). 

Proposition 5.2 illustrates that we may indeed incorporate deCastel- 
jau's algorithm when evaluating mp(X s ). 
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PROPOSITION 5.2. Let n > s > 1, (3 G Z s + , and X = {x°, . . . , x n } C R s + 
with vol s ([X}) > 0. Let y i = (a£, . . . ,3?) G K+ +1 and set g l = \y% 
i = 1, . . . , s. Then 

to *W = MtU ^ E (5.15) 

J ' |fci|=0i |fc.|=0. 

s 

Here, h = (&»„,. . . , k in ) G Z+ +1 , 77! = r] l ■ ■ -r] n l, with rjj = Yl hj, 

i=i 

0<j<n,y l = (l/g l )y\ i = 1, . . . , s. 
9 = (go,---,9n) 

Proof. We use the defining relation (2.2) to write: 
mp(X) =n\ I f\(t x° 3 +■■■ + t n x])h dt 

= *, E (*)(✓>»■•■• E gW 



where 77, fcj, i = 1, . . . , s, are given in Proposition 5.2. Using 4.3-4 in [6] 
to simplify the integral in the above identity, we have 

Now y J , j = 1, . . . , s, given in Proposition 5.2 can be viewed as the 
barycentric coordinates of some point in M. n . Inserting y J, s into (5.16) and 
scaling by g 13 gives the desired result. □ 

Thus all moments of the form mp{X s ) can be expressed as a nested 
sum of Bezier polynomials and subsequently may be evaluated using de- 
Casteljau's algorithm. It should be noted that while both de Casteljau's algo- 
rithm and (5.14) are possible candidates for the task of computing m^X^), 
k = s + 1, . . . , n, the latter scheme requires the evaluation of fewer terms at 



each recursion step and is thus the preferred choice. In order to summarize 
the procedure for evaluating the moments of simplex splines in Algorithm 
5.3, we introduce d = (0, . . . , 0) G R s . 

ALGORITHM 5.3. 

Given X = {x ,... ,x n } C R s , (n > s) and f3 e Z s + with \/3\ > 1, this 

algorithm generates the moment mp{X) of the simplex spline M(-\X). 

1. a := do 

2. For k = s to n 

I — ► m a (X k ) = 1 

3. For \a\ = 1 to \/3\, aeZ s + ,a<f3 

I > Use (5.15) to express m a (X s ) in terms of Bezier polynomials 

and evaluate using DeCasteljau's algorithm. 

4. For k = s + 1 to n 

| For a e Z s + , a < (3 

I >. I > Compute m a (Xk) using (5.14). 

We close this section with a remark that this algorithm is numerically 
stable if x J > for all j = 0, 1, . . . , n. 

6. Applications to Hypergeometric Functions. 

In this section we demonstrate a relationship between Dirichlet splines 
and an important class of hypergeometric functions of several variables. We 
will deal mainly with Appell's F 4 and Lauricella's F B . The link between 
these classes of functions is provided by another integral average which is 
commonly referred to as a double Dirichlet average (see [4] for more de- 
tails). Throughout the sequel the double Dirichlet average of a continuous 
univariate function h will be denoted by 7i. 

Let X e M sx ( n+1 )(n > s > 1). Further, let u = (u u . . . ,u 8 ) be an 
ordered s-tuple of nonnegative numbers with u\ + • • • + u s = 1, and similarly 



v = (vo, . . . , v n ). We define 

s n 

It • Xv = UiX^Vj, 
i=l j=0 

where x^ stands for the i-th component of the j-th column of X . Let h be a 
continuous function on i" = [Min x\, Max x?]. In order to avoid trivialities, 
we will assume that I has a nonempty interior. For b = (pi, . . . , b 8 ) e 
and d = (d , • • • , d n ) E M> +1 , let [4, p. 421] 



H(b;X;d) = 




h(u ■ Xv)<f>b(u)<f)d(v) dudv, 



du = du2 ■ ■ ■ du s , dv = dv\ ■ ■ ■ dv n . Here <f>b and <f>d are the Dirichlet densities 
on E s _i and E n , respectively (see (2.1)). It is known that for b G 

H(b;X;d)= I H(d;u-X)(j) b (u)du 

(see [4, (2.8)]). In (6.1) H stands for the single Dirichlet average of h (see 
(3.1)), u • X = {u ■ x°, . . . , u ■ x n }, x°, . . . , x n - the columns of X. 

We are in a position to state and prove the following: 

THEOREM 6.1. Let d € M> +1 and let the vector b G W be such that 
c^O, —1, ... ,(c = bi-\ + b a ). If vol 8 ([X}) > 0, then 

H(b;X;d)= I M(x\d;X)H(b;x)dx, 

J[X] 

x = (xi, . . . , x s ), dx = dxi ■ ■ ■ dx s . 

Proof. In order to establish (6.2) assume for the moment that b G M>. Ap- 
plication of (3.2) and (3.10) to (6.1) gives 

H(b;X;d)= [ h(u ■ x)M(x\d; X) dx]<f> b (u) du. 

iE a _! J[X] 
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Interchanging the order of integration and next using (3.1), we obtain 
the assertion provided b G R>. This restriction can be dropped because the 
average H can be continued analytically in the 6-parameters provided that 
c 7^ 0, —1, ... , (see [6,Thm. 6.3-7]). This completes the proof. □ 

Before we state a corollary of Theorem 6.1 let us introduce more 
notation. For h(z) = z~ a , (a G R), the double Dirichlet average of h will be 
denoted by TZ- a (cf. [4]). 

COROLLARY 6.2. ([7]) Let d G R> +1 , b G R s , and let the matrix X be such 
that S £ [X]. Then 

m- b (d;X) = K- c (b;X;d), (6.3) 

where m-b(d; X) stands for the moment of order — c of the Dirichlet spline 
M(-\d;X). 

Proof. Apply [6,(6.6-5)] 

s 

R- c (b;X) = l[x- b * (6.4) 

i=l 

to (6.2) with h(t) = t~ c . □ 

Hereafter, we will deal with the hypergeo metric functions and poly- 
nomials of several variables. Appell's hypergeometric function F 4 is defined 
by the double power series [6, Ex. 6.3-5] 

F 4 (a 1 P;-f,5;x 1 ,X2) = > > — — .., x . v , ., — x x x 2 \ 

7~^f~i (7, WW 

a,/?,7,(5 G R, 7, 6 ^ 0,-1,..., \xi\? + \x2\^ < 1. The following integral 
formula [5, p. 963] 

F 4 (a,/9;7,<y;xi(l -x 2 ),x 2 (l -xi)) = 



f 

Jo 



R- a (do, di, d 2 ] u ■ x , -u • x 1 , -u • rr 2 )^^) du (6.5) 
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provides the analytic continuation of the F 4 -series to the region A defined 
by 

A = {(xi, x 2 ) G M 2 : x\ < 1, x 2 < 1, xi + x 2 < 1}. 

In (6.5), b = (/?, 7 — /?), d = j + 6-a-l, di = a + fi-^-S+l, d 2 = 6-0, 
4>b{u) is the Dirichlet density on E\, and x° , x 1 , and x 2 are the columns of 
X, where 

(1 - Xi){l — x 2 ) \ — x\—x 2 1—X\ 
1 - x 2 1 - x 2 1 



X = 



(6.6) 



COROLLARY 6.3. Let d = (d , d u d 2 ) G R> and let b = (p, 7 - (3) G K 2 . If 
vol 2 {[X\) > 0, then 

F 4 (a,fc 1 ,8;x 1 (l-x 2 ),x 2 (l-x 1 )) = [ M(y\d; X)R_ a (b;y) dy (6.7) 

J[X] 

(y = (yi,y 2 ),dy = dy\dy 2 ). Here, R- a is the single Dirichlet average of 
h(z) = z~ a and the matrix X is given in (6.6). 

Proof. Apply (6.1) to (6.5) and next use (6.2). □ 
A special case of (6.7) is 

F 4 (o!,/?; a, 8; xi(l —x 2 ),x 2 (l - xi)) = m-b(d;X), 

where now b = (/?, a — (3) and d = (8 — 1, (3 — 8 + 1, 8 — (3). This follows 
immediately from (6.7) and (6.4). 

We will now deal with Lauricella's Fb function and Lauricella poly- 
nomials. Let a = (ai, . . . , a n ) G R n , (3 = (fa, . . . , (3 n ) G R n , 7 G R 
(7 7^ 0,-1,...), and let x = (xi,...x n ) G R n , with \xi\ < 1, for all i. 
Following [21] we define 

F B {a, (3; 7; x) = ^ ( 7? | fc |) fc! x ' ( 6 - 8 ) 
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where the summation extends over all multi-indices k = (k±, . . .k n ) £ Z™. 
In (6.8) we employ multi-index notation introduced in Section 2. Also, 

n 

(a,k) = JJ(«i,/ci). 

i=i 

(/3, fc) is defined in an analogous manner. When n = 1, becomes Gauss' 
2-F1 function. 

COROLLARY 6.4. Let d = 7 - e M> +1 and let 



X 



1-xi 1 

1 1-X 2 



1 1 
1 1 



Then 



1 ... 1 - x n 1 



F B (a,(3-^]x) = m- a (d;X) 



(6.9) 



(6.10) 



provided that Xi < 1 for all i. 

Proof. In the stated domain the entries of X are positive. Thus n ^ [X]. 
In order to establish (6.10) we utilize [4,(5.11)] to obtain 

F B (a;(3;-f;x) = 1Z_ c (a;X;d), 

where now c = ol\ + • • • + ct n . This in conjunction with (6.3) gives the 
assertion. □ 

Lauricella polynomials Lj(x), (j G Z^,x e M n ) are defined in the 
following way: 

These polynomials play an important role in the study of coherent states (cf. 
[16]). On account of (6.10) 



Lj{x) = rrij(d; X) 



(6.11) 
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where the vector d and the matrix X are the same as in Corollary 6.4. 

Two generating functions for the polynomials under discussion can 
be derived from (5.9) and (5.11). Let A = (Ai,... , A n ), e = (1,... ,1), 
A,e e M n , and let d= (/?,7 - \(3\) E R n+1 . Then 

exp(X ■ e)S(d; -Xixi, . . . , -A n ;r n ,0) = ^ —Lj(x) (6.12) 

j e z\. if 

Max{|A • e — Aixi|, ... , |A • e — A n x n |, |A • e\} < 1, 

then 

R- a (d-Y) = Y^X^p-L j {x\ (6.13) 

a e R, j e Z^, Y = {1 - A • e + Aixi, ... , 1 - A • e + A n x n , 1 - A • e}. 

For the proof of (6.12) we replace b by d in (5.9) and next use (6.11), 
(6.9), and [6,(5.8-3)]. This gives the desired result provided that d G R> +1 . 
The latter restriction can be dropped because the S- function can be con- 
tinued analytically in the <i-parameters ([6, Corollary 6.3-3]). The generating 
function (6.13) can be derived from (5.11) by the same means. Feinsilver's 
generating function [15] can be obtained from (6.13) by letting a = 7 and 
then using (4.5). Recall that R = 1. It is not hard to show that the R- 
hypergeometric function R- a in (6.13) is a multiple of the Lauricella function 
of the fourth kind. We have 



R_ a (d; Y) = (1 - A • e)- a F D (a, (3; 7; z h . . . , z n ), 

where 



AjXj 

Zi ~ - 



1- A-e' 



i = 1, 2, . . . , n. We omit further details. 
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Our next goal is to establish a recurrence formula obeyed by 
Lauricella polynomials 

(7 + \k\)L k+dm (x) - [7(1 - w m x m ) + \k\]L k (x) 

n 

+ ^2 kl6lm \- Lk ~ d i ( x ) ~~ L k-d l +d m {x)\ = (6.14) 
1=1 

m = 1, 2, . . . ,n. In (6.14), 7 G R>, fe = (fci, . . . , fc n ) G Z™, <i m stands for the 
mth coordinate vector in IR n , similarly di, w m = (3 m > 0, 1 < m < n, 

fim= 1 ., . (6-15) 

[ 1 — x m 11 I = m. 

Here we adopt the convention that L k (x) = if — k m G N for some m. In 
order to establish the recursion (6.14) we derive first a recurrence formula 
for the moments of multivariate Dirichlet splines with s = n and X given in 
(6.9). We have 

(c + \p\) mi3+dm (6; X) - [c(l - w m x m ) + \[3\}mp(b; X) 

n 

+ ^2(3ie lm [mp_ dl (b;X) - mp- dl+dm (b; X)] = 0, (6.16) 
1=1 

beR n+1 ,c = b! + --- + b n+1 , P eR n , \/3\ = ft + •••/?„. The recursion (6.14) 
now follows from (6.16) by letting b = d, f3 = k G Z™, and using (6.11). 
To complete the proof we need to establish (6.16). To this aim we increase 
the indices of summation in (5.2) and (5.3) by one unit. Next we let s = n 
and solve the resulting linear system for mp(b + e m ;X), 1 < m < n + 1. 
Let us note that the assumption vol n ([X]) > is equivalent to x m 7^ 0, 
1 < m < n. This assures uniqueness of the solution. Subtracting (5.3) from 
(5.2) we obtain 

m p (b + e m - X) = ^p(b;X)-m, +dm (b;X)] ^ ^ 
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1 < m < n. The remaining moment mp(b + e n+ i]X) can be found using 
(5.2) and (6.17). To complete the proof of (6.16) we utilize (5.4). Replacing 
the index j by m and next using (6.17), we can easily obtain the assertion. 

We close this section with an inequality for Lauricella polynomials. 
To this end, let x t < 1, 1 < i < n. It follows from (6.9) that [X] C R> in 
the stated domain. This in conjunction with (6.11) and (5.1) provides 

Lj(x) = [ yJM(y\d;X)dy>0 

J[X] 

(j ez»,d= e k> +1 , y = (yu--- dy = dy u ...dy n ). A 

standard argument applied to the last formula gives 

[Lj(x)] 2 < Lj_ k {x)L j+k {x), 

where the vector k G is such that j — k e Z". In particular, if 
k = e m - the mth coordinate vector in M n , then 

[L 3 (x)f < Lj- em (x)L j+em (x), 

provided j — e m G Z" . Thus the function g : Z" — > R, where 

is log-convex in each variable separately. 
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"Moments of Dirichlet Splines and Their 
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by 
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K. upper case open-face Roman R 

> greater than or equal to 
G is a member of 

summation sign 

C is a subset of 

N upper case open-face Roman N 

Z upper case open-face Roman Z 

(3 lower case Greek beta 

! factorial 

a lower case Greek alpha 

< less than or equal to 

> greater than 

<j> lower case Greek phi 

f | product 

T upper case Greek gamma 

f integral sign 

T upper case script F 

H upper case script H 

1Z upper case script R 

e lower case Greek epsilon 

A lower case Greek lambda 

7 lower case Greek gamma 

5 lower case Greek delta 

partial derivative 

A upper case Greek lambda 
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